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Abstract 

By measuring hadronic single-particle spectra and two-particle correla¬ 
tions in heavy-ion collisions, the size and dynamical state of the collision 
fireball at freeze-out can be reconstructed. I discuss the relevant the¬ 
oretical methods and their limitations. By applying the formalism to 
recent pion correlation data from Pb+Pb collisions at CERN we demon¬ 
strate that the collision zone has undergone strong transverse growth 
before freeze-out (by a factor 2-3 in each direction), and that it expands 
both longitudinally and transversally. From the thermal and flow en¬ 
ergy density at freeze-out the energy density at the onset of transverse 
expansion can be estimated from conservation laws. It comfortably ex¬ 
ceeds the critical value for the transition to color deconfined matter. 


1 Introduction 

In the last few years a large body of evidence has been accumulated that the 
hot and dense collision region in ultrarelativistic heavy ion collisions thermal- 
izes and shows collective dynamical behaviour. This evidence is based on a 
comprehensive analysis of the hadronic single particle spectra. It was shown 
that all available data on hadron production in heavy ion collisions at the AGS 
and the SPS can be understood within a simple model which assumes locally 
thermalized momentum distributions at freeze-out, superimposed by collective 
hydrodynamical expansion in both the longitudinal and transverse directions 
[;1], [2j. The collective dynamical behaviour in the transverse direction is re¬ 
flected by a characteristic dependence of the inverse slope parameters of the 
mj_-spectra (“effective temperatures”) at small m± on the hadron masses [|lj]. 
New data from the Au+Au and Pb+Pb systems || support this picture and 
show that the transverse collective dynamics is much more strongly exhibited 
in larger collision systems than in the smaller ones from the first rounds of 
experiments. The amount of transverse flow also appears to increase mono- 
tonically with collision energy from GSI/SIS to AGS energies, but may show 


signs of saturation at the even higher SPS energy 0. 

The extraction of flow velocities and thermal freeze-out temperatures from 
the measured single particle spectra relies heavily on model assumptions [[1]. 
There have been alternative suggestions to explain the observed features of the 
hadron spectra without invoking hydro-flow 0, |5j . The single-particle spectra 
are ambiguous because they contain no direct information on the space-time 
structure and the space-momentum correlations induced by collective flow. In 
terms of the phase-space density at freeze-out (“emission function”) S(x,p ) 
the single-particle spectrum is given by EdN/d 3 p = f d 4 x S(x, p); the space- 
time information in S is completely washed out by integration. Thus, on 
the single-particle level, comprehensive model studies are required to show 
that a simple hydro dynamical model with only a few thermodynamic and 
collective parameters can fit all the data, and additional consistency checks are 
needed to show that the extracted fit parameter values lead to an internally 
consistent theoretical picture. The published literature abounds with examples 
demonstrating that without such consistency checks the theoretical ambiguity 
of the single particle spectra is nearly infinite. 

This is the point where Bose-Einstein correlations between the momenta 
of identical particle pairs provide crucial new input. They give direct access to 
the space-time structure of the source and its collective dynamics. In spite of 
some remaining model dependence the set of possible model sources can thus be 
reduced dramatically. The two-particle correlation function C(q,K) is usually 
well approximated by a Gaussian in the relative momentum q whose width 
parameters are called “HBT (Hanbury Brown-Twiss) radii”. It was recently 
shown P]-[Q that these radius parameters measure certain combinations of 
the second central space-time moments of the source. In general they mix the 
spatial and temporal structure of the source in a nontrivial way [ij7|], and the 
remaining model dependence enters when trying to unfold these aspects. 

Collective dynamics of the source leads to a dependence of the HBT radii 
on the pair momentum K: this has been known for many years [TH, ITT], but 


was recently quantitatively reanalyzed, both analytically 1B i |E| and nu¬ 
merically |L3], [L4|. The velocity gradients associated with collective expansion 
lead to a dynamical decoupling of different source regions in the correlation 
function, and the HBT radii measure the size of the resulting “space-time 
regions of homogeneity” of the source 00 around the point of maximum 
emissivity for particles with the measured momentum K. The velocity gradi¬ 
ents are smeared out by a thermal smearing factor arising from the random 
motion of the emitters around the fluid velocity dj. Due to the exponential 
decrease of the Maxwcll-Boltzmann distribution, this smearing factor shrinks 













with increasing transverse momentum K± of the pair, which is the basic reason 
for the /ix-dependence of the HBT radii. 

Unfortunately, other gradients in the source (for example spatial and tem¬ 
poral temperature gradients) can also generate a /t-dependence of the HBT 
radii 0, [L2|. Furthermore, the pion spectra in particular are affected by res¬ 
onance decay contributions, but only at small Jij_- This may also affect the 
HBT radii in a AA-dependent way [[To|. [H| . The isolation of collective flow, 


in particular transverse flow, from the AA-dependence of the HBT radii thus 
requires a careful study of these different effects. 

We studied this AT-dependence of the HBT radii within a simple analytical 
model for a finite thermalized source which expands both longitudinally and 
transversally. For presentation 1 use the Yano-Koonin-Podgoretskii (YKP) 
parametrization of the correlator which, for sources with dominant longitu¬ 
dinal expansion, provides an optimal separation of the spatial and temporal 
aspects of the source || [T4| |. The YKP radius parameters are independent of 
the longitudinal velocity of the observer frame. Furthermore, in all thermal 
models without transverse collective flow, they show perfect M ±-scaling (in 
the absence of resonance decay contributions). Only the transverse gradients 
induced by a non-zero transverse flow can break this Mj_-scaling, causing an 
explicit dependence on the particle rest mass. This allows for a rather model- 
independent identification of transverse flow from accurate measurements of 
the YKP correlation radii for pions and kaons. High-quality data should also 
allow to control the effects from resonance decays. 

Due to space limitations, 1 will be extremely selective with equations, fig¬ 
ures and references. A comprehensive and didactical discussion of the for¬ 
malism and a more extensive selection of numerical examples can be found in 
Ref. 


17 to which I refer the reader for more details. 


2 The source model 

For our study we use the model of Ref. [§] for an expanding thermalized source: 

K-u{x) (r - r 0 ) 2 r 2 (77 - r/ 0 ) 2 
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( 1 ) 


Here r 2 = x 2 + y 2 , the spacetime rapidity y = | ln[(t + z)/(t — z)], and the 
longitudinal proper time r = \Jt 2 — z 2 parametrize the spacetime coordinates 
oY, with measure d A x = t dr dyr dr d(j). Y — | In [(Ek + Kl )/ [Ek — Kl)\ and 

M± = Jm 2 -\-K\ parametrize the longitudinal and transverse components 



















of the pair momentum K. T(x) is the freeze-out temperature, a/2 R is the 
transverse geometric (Gaussian) radius of the source, r 0 its average freeze-out 
proper time, At the mean proper time duration of particle emission, and Aq 
parametrizes 0 the finite longitudinal extension of the source. The expansion 
flow velocity u^(x) is parametrized as 

u^(;x) = (cosh?/cosh rj t (r), sinh rj t (r) e r , sinh?/cosh?//r)), (2) 

with a boost-invariant longitudinal flow rapidity rji = rj (vi = z/t) and a linear 
transverse flow rapidity profile 

Vt(r) = T) f (-0 . (3) 

rjf scales the strength of the transverse flow. The exponent of the Boltzmann 
factor in ([!]) can then be written as 

K ■ u(x ) = M± cosh(y — rj) cosh rj t {r) — K±-e r sinh rj t (r) . (4) 

For vanishing transverse flow (rjf = 0) the source depends only on M±, and 
remains azimuthally symmetric for all K±. 


3 HBT radius parameters 


From the source function (|T|) the correlation function is calculated via the 


relation [IS 


C(q, K) « 1 + 


| Jd 4 xS(x,K) e iq x 
[/ d 4 x S(x, K)\ 2 


= 1 + 


fd 4 xS(x,K) eW-fa 


\Jd 4 xS(x,K)\' 


( 5 ) 


Here q = P\—p 2 ) K = (P 1 +P 2 )/2 are the relative and average 4-momenta of the 
boson pair. The quality of the approximation in (|]) is discussed in 0. Since 
Pi,P 2 are on-shell and thus K ■ q — 0, q° = (3-q (with f3 = K / K° ~ K/Ek) 
is not an independent variable and can be eliminated (second equality in (|5D). 
Therefore the Fourier transform in (|5|) cannot be inverted without a model for 
S(x,K). This is the reason for the model dependence of the interpretation of 
HBT correlation data mentioned in the Introduction. 

We use a cartesian coordinate system with the z-axis along the beam di¬ 
rection and the x-axis along K±, Then (3 = (/3±,0,/3i). We assume an az¬ 
imuthally symmetric source (impact parameter ss 0) and write C(q,K) in the 
YKP parametrization |T4||: 


C(q, K) 


1 + exp 
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( 6 ) 











Here q'j_ — q 2 + q 2 , and R±, Ay, Ro, U are four A'-dependent parameter func¬ 
tions. U(K ) is a 4-velocity with only a longitudinal spatial component: 


U(K) = 7 (K) (l, 0, 0, v(K)) , with 7 = . (7) 

Its value depends, of course, on the measurement frame. The “Yano-Koonin 
velocity” v(K) can be calculated in an arbitrary reference frame from the 
second central space-time moments of S(x,K ) (for explicit expressions see 
Ref. [|hj). It is, to a good approximation, the longitudinal velocity of the 
fluid element from which most of the particles with momentum K are emitted 
[§, [b|. For sources with boost-invariant longitudinal expansion velocity the 
YK-rapidity associated with v(K) is linearly related to the pair rapidity Y 
0 - 

The other three YKP parameters do not depend on the longitudinal velocity 
of the observer. (This distinguishes the YKP form (^j) from other parametriza- 
tions 0 !H0 Their physical interpretation is easiest in terms of coordinates 
measured in the frame where v(K) vanishes. There they are given by j8| 


R\(K) = 

(: v 2 >, 

(8) 

R\(K) = 

((5 - (f3i/f3 ± )x)' 2 ) - ( (3i/(3±) 2 (y 2 ) « (z 2 ), 

(9) 

Rl(K) = 

/(t-x/pA 2 ) - (y 2 )/(3 2 ± w ( i 2 ). 

(10) 


Here ( f(x )) = / d 4 x f(x) S(x, K)j J d 4 x S(x, K) denotes the (A"-dependent) 
average over the source function S(x,K), and x = x — x(K) etc., where 
x{K) = ( x ) is (approximately) the source point with the highest intensity 
at momentum K. Aj_, R» and Ro thus measure, approximately, the (In¬ 
dependent) transverse, longitudinal and temporal regions of homogeneity of 
the source in the local comoving frame of the emitter. The approximation in 
(IDi JTOl) consists of dropping terms which (for our model) vanish in the absence 
of transverse flow and were found in |14| to be small even for finite transverse 
flow. Note that it leads to a complete separation of the spatial and temporal 
aspects of the source. This separation is spoiled by sources with (x 2 ) ^ (x 2 ). 
For our source this happens for non-zero transverse (in particular for large) 
transverse flow rjf, but for opaque sources where particle emission is surface 


dominated [|l9j this may be true even without transverse flow. 

Since in the absence of transverse flow the /^-dependent terms in (I9|) and 
© vanish and the source itself depends only on M±, all three YKP radius 




















parameters then show perfect Mj_-scaling. Plotted as functions of M±, they 
coincide for pion and kaon pairs (see Fig. 5 in |L7|). For non-zero transverse 
flow this Mj_-scaling is broken by two effects: (1) The thermal exponent (£f|) 
receives an additional contribution proportional to K± = JM\ — m 2 . (2) The 
terms which were neglected in the second equalities of ([^.[fl]) are non-zero, and 
they also depend on f3± = K^/E k . Both effects induce an explicit rest mass 
dependence and destroy the Mj_-scaling of the YKP size parameters. 


4 My-dependence of YKP radii and collective flow 


Collective expansion induces corelations between coordinates and momenta in 
the source, and these result in a dependence of the HBT parameters on the 
pair momentum K. At each point in the source the local velocity distribution 
is centered around the average fluid velocity; two points whose fluid elements 
move rapidly relative to each other are thus unlikely to contribute particles with 
small relative momenta. Essentially only such regions in the source contribute 
to the correlation function whose fluid elements move with velocities close to 
the velocity of the observed particle pair. If the source expands rapidly and 
features large velocity gradients, these contributing “regions of homogeneity” 
will be small. Their size will be inversely related to the velocity gradients, 
scaled by a “thermal smearing factor” ^Jt/M± which characterizes the width 
of the Boltzmann distribution [!|7]. 

Thus, for expanding sources, the HBT radius parameters are generically 
decreasing functions of the transverse pair mass M±. The slope Of this de¬ 
crease grows with the expansion rate [|1^, [Uj. Longitudinal expansion affects 
mostly the longitudinal radius parameter R» and the temporal parameter R 0 
|Tl|j ; the latter is a secondary effect since particles from different points are 
usually emitted at different times, and a decreasing longitudinal homogene¬ 
ity length thus also leads to a reduced effective duration of particle emission. 
The transverse radius parameter R± is invariant under longitudinal boosts and 
thus not affected at all by longitudinal expansion. It begins to drop as a func¬ 
tion of M_l, however, if the source expands in the transverse directions. The 
sensitivity of and Rq to transverse flow is much weaker [Q. Transverse 
(longitudinal) flow thus affects mostly the transverse (longitudinal) regions of 
homogeneity. For a quantitative study see Fig. 6 in Ref. 0- 

Unfortunately, the observation of an Mj_-dependence of R± by itself is 
not sufficient to prove the existence of radial transverse flow. It can also be 
created by other types of transverse gradients, e.g. a transverse temperature 
gradient [0, O. To exclude such a possibility one must check the Mj_-scaling 






of the YKP radii, i.e. the independence of the functions Ri(M± ) ( i =_L, ||, 0) of 
the particle rest mass (which is not broken by temperature gradients). Since 
different particle species are affected differently by resonance decays, such a 
check further requires the elimination of resonance effects. 


5 Resonance decays 


Resonance decays contribute additional pions at low Mj_; these pions originate 
from a larger region than the direct ones, due to resonance propagation before 
decay. They cause an Mj_-dependent modification of the HBT radii. 

fl have shown that the resonances can be sub di- 


Quant it at ive studies |15 


vided into three classes with different characteristic effects on the correlator: 

(i) Short-lived resonances with lifetimes up to a few fm/c do not propagate 
far outside the region of thermal emission and thus affect R± only marginally. 
They contribute to Ro and R\\ up to about 1 fin via their lifetime; R\\ is larger 
if pion emission occurs later because for approximately boost-invariant expan¬ 
sion the longitudinal velocity gradient decreases as a function of time. 

(ii) Long-lived resonances with lifetimes of more than several hundred fm/c 
do not contribute to the measurable correlation and thus only reduce the cor¬ 
relation strength (the intercept at q — 0), without changing the shape of the 
corelator. The reason is that they propagate very far before decaying, thus 
simulating a very large source which contributes to the correlation signal only 
for immeasurably small relative momenta. 

(iii) There is only one resonance which does not fall in either of these two 
classes and can thus distort the form of the correlation function: the u with 
its lifetime of 23.4 fm/c. It contributes a second bump at small q to the corre¬ 
lator, giving it a non-Gaussian shape and complicating Jl6| the extraction of 
HBT radii by a Gaussian fit to the correlation function. At small Mi up to 
10% of the pions can come from a; decays, and this fraction doubles effectively 
in the correlator since the other pion can be a direct one; thus the effect is not 
always negligible. 

In a detailed model study |L6| we showed that resonance contributions can 
be identified through the non-Gaussian features in the correlator induced by 
the tails in the emission function resulting from resonance decays. To this 
end one computes the second and fourth order g-moments of the correlator 
[llffl . The second order moments define the HBT radii, while the kurtosis 
(the normalized fourth order moments) provide a lowest order measure for 
the deviations from a Gaussian shape. We found (16] that, at least for the 


model ([[]), a positive kurtosis can always be associated with resonance decay 




















contributions. Strong flow also generates a non-zero, but small and apparently 
always negative kurtosis (see Fig. 12 in 0]). Any M^-dependence of R± which 


is associated with a positive Mj_-dependent kurtosis must therefore be regarded 
with suspicion; an M ±-dependence of R± with a vanishing or negative kurtosis, 
however, cannot be blamed on resonance decays. 



Fig. 1. The inverted (/-variance R± and the kurtosis A_l (the index s in the 
figure stands for “sideward”) at Y = 0 as functions of K±_. Left: r/j = 0 (no 
transverse flow). Right: r)f = 0.3. The difference between dashed and solid 
lines is entirely dominated by u decays. 

In our model, the first situation is realized for a source without transverse 
expansion (left panel of Fig. 1): At small M± the oj contribution increases 
R _l by up to 0.5 fin while for Mj_ > 600 MeV it dies out. The effect on R± 
is small because the heavy u moves slowly and doesn’t travel very far before 
decaying. The resonance contribution is clearly visible in the positive kurtosis 
(lower curve). For non-zero transverse flow (right panel) there is no resonance 
contribution to R±] this is because for finite flow the effective source size for the 
heavier c o is smaller than for the direct pions, and the cn-decay pions thus always 
remain buried under the much more abundant direct ones. Correspondingly 
the kurtosis essentially vanishes; in fact, it is slightly negative, due to the weak 
non-Gaussian features induced by the transverse flow. 

6 Analysis of Pb+Pb data 

In Fig. 2 we show a numerical fit of the YKP radius parameters, using the 
expressions (D)-(|l(]|) with our model source ([I]), to data collected by the NA49 
















collaboration in 158 A GeV/c Pb+Pb collisions [^, |2l| . Please note that this 
fit refers to only a single rapidity slice of the available data, and it does not 
include resonance decays. The fit result must therefore be taken with great 
care. A comprehensive simultaneous analysis of all single particle spectra and 
two-particle correlation data from Pb+Pb collisions is in progress. 


Data: Pb+Pb 158 AGeV, 
NA49 preliminary, 
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Fig. 2. R l, R\\ and Rq for 158 A GeV/c Pb+Pb collisions as functions of 
the transverse pair momentum. The data are from the NA49 Collaboration 
|2lj |. The lines are a fit with the model (ifj), with fit parameters as given in the 
figure. 


After A?7=1.2 has been adjusted to reproduce the width of the pion rapidity 
distribution |j2lJ , the parameters To and Ar] are essentially fixed by the magni¬ 
tude of R\\ and i?o- The radius R is fixed by the magnitude of R±(K± = 0) once 
the temperature T and transverse flow rjf are known. The Mj_-dependence of 
















R l fixes T and rjf , albeit not independently: essentially only the combination 
r]f \JM±/T , i.e. the velocity gradient divided by the thermal smearing factor, 
can be extracted P, [2lJ. This is similar to the single particle spectra whose 


M ±-slopes determine only an effective blushifted temperature, T e g = T\j 

fl[]. The correlations between T and rjf are, however, exactly opposite in the 
two cases: for a fixed spectral slope T must be decreased if rjf increases while a 
fixed M_L-slope of R± requires decreasing values of rjf if T is reduced [|2T]]. The 
combination of single-particle spectra and two-particle correlation thus allows 
for a separate determination of T and rjf. 

For the fit in Fig. 2 the freeze-out temperature was set by hand to T = 
100 MeV. The resulting flow parameter rjf— 0.6 corresponds to an average 
transverse flow velocity b/=0.58. This combination of T and rjf results in 
single-particle spectra with roughly the right shape. Fitting R± with higher 
temperatures results in larger rjf- values which leads to single particle spectra 
which are much too flat. 

Let us discuss in more detail the numbers resulting from this fit. First, the 
transverse size parameter R—7 fm is surprisingly large. Resonance contribu¬ 
tions are not expected to reduce it by more than 0.5 fm |TB[. The transverse 


flow correction to R± is appreciable, resulting in a visible transverse homogene¬ 
ity length of only about 5.5 fm at small K±, but even this number is large. 
R—7 fm corresponds to an r.m.s. radius r rms = J ( x 2 + y 2 ) ~ 10 fm of the 


pion source, to be compared with an r.m.s. radius r^ s = 1.2 x TL 1 / 3 /-^ fm = 
3.2 fm for the density distribution of the original Pb nucleus projected on the 
transverse plane. This implies a transverse expansion of the reaction zone by 
a linear factor 3. That we also find a large transverse flow velocity renders the 
picture consistent. The longitudinal size of the collision region at the point 
where the pressure in the system began to drive the transverse expansion is less 
clear; assuming a similar expansion factor 3 in the beam direction (although 
we know from the fit that the longitudinal expansion velocity is much bigger 
than the transverse one) we conclude that the fireball volume must have ex¬ 
panded by at least a factor 3 3 = 27 between the onset of transverse expansion 
and freeze-out! This is the clearest evidence for strong collective dynamical 
behaviour in ultra-relativistic heavy-ion collisions so far. 

The local comoving energy density at freeze-out can be estimated from the 
fitted values for T and rjf. The thermal energy density of a hadron resonance 
gas at T = 100 MeV and moderate baryon chemical potential is of the order 
of 50 MeV/fm 3 . The large average transverse flow velocity of (vf) ~ 0.58 
implies that about 50% flow energy must be added in the lab frame. This 















results in an estimate of about 0.050 GeV/fm 3 x 1.5 x 27 ~ 2 GeV/fm 3 for the 
energy density of the reaction zone at the onset of transverse expansion. This 
is well above the critical energy density e cr < 0.9 GeV/fm 3 predicted by lattice 
QCD for deconhned quark-gluon matter |22|]. Whether this energy density was 
thermalized is, of course, a different question. In any case it must have been 
accompanied by transverse pressure, because otherwise transverse expansion 
could not have been initiated. 

7 Conclusions 

I hope to have shown that 

• two-particle correlation functions from heavy-ion collisions provide valu¬ 
able information both on the geometry and the dynamical state of the 
reaction zone at freeze-out; 

• a comprehensive and simultaneous analysis of single-particle spectra and 
two-particlc correlations, with the help of models which provide a real¬ 
istic parametrization of the emission function, allows for an essentially 
complete reconstruction of the final state of the reaction zone, which can 
serve as a reliable basis for theoretical back-extrapolations towards the 
interesting hot and dense early stages of the collision; 

• simple and conservative estimates, based on the crucial new information 
from HBT measurements on the large transverse size of the source at 
freeze-out and using only energy conservation, lead to the conclusion 
that in Pb+Pb collisions at CERN, before the onset of transverse expan¬ 
sion, the energy density exceeded comfortably the critical value for the 
formation of a color deconhned state of quarks and gluons. 

Acknowledgements: This work was supported by grants from DAAD, DFG, NSFC, 
BMBF and GSI. The results reported here were obtained in collaboration with S. Chap¬ 
man, P. Scotto, B. Tomasik, U. Wiedemann, and Y.-F. Wu, to whom I would like express 
my thanks. I gratefully acknowledge many discussions with H. Appelshauser, T. Csorgo, D. 
Ferenc, M. Gazdzicki, S. Schonfelder, P. Seyboth, and A. Vischer. 

References 

[1] K.S. Lee and U. Heinz, Z. Phys. C43 (1989) 425; K.S. Lee, U. Heinz and 
E. Schnedermann, Z. Phys. C48 (1990) 525; E. Schnedermann and U. 
Heinz, Phys. Rev. Lett. 69 (1992) 2908; E. Schnedermann, J. Sollfrank 
and U. Heinz, Phys. Rev. C48 (1993) 2462; E. Schnedermann and U. 


Heinz, Phys. Rev. C50 (1994) 1675; U. Heinz, in Hot Hadronic Matter: 
Theory and Experiment , Eds: J. Letessier et al., NATO ASI Series B346, 
Plenum, New York, 1995, p.413. 

[2] J. Stachel et al., [E814 collaboration], Nucl. Phys. A566 (1994) 183c; P. 
Braun-Munzinger et al., Phys. Lett. B344 (1995) 43. 

[3] See proceedings of Quark Matter 96, Nucl. Phys. A610 (1996), in press (in 
particular contributions by Y. Akiba, N. Herrmann, P. Jones, R. Lacasse, 
J. Stachel, and Nu Xu). 

[4] A. Leonidov, M. Nardi and H. Satz, Z. Phys. C, in press. 

[5] S. Esumi, U. Heinz, and N. Xu, submitted to Phys. Lett. B. 

[6] M. Herrmann and G.F. Bertsch, Phys. Rev. C51 (1995) 328. 

[7] S. Chapman, P. Scotto and U. Heinz, Phys. Rev. Lett. 74 (1995) 4400; 
and Heavy Ion Physics 1 (1995) 1. 

[8] S. Chapman, J.R. Nix and U. Heinz, Phys. Rev. C52 (1995) 2694. 

[9] S.V. Akkelin and Y.M. Sinyukov, Phys. Lett. B356 (1995) 525. 

[10] S. Pratt, Phys. Rev. Lett. 53 (1984) 1219; Phys. Rev. D33 (1986) 1314. 

[11] A.N. Makhlin and Yu.M. Sinyukov, Z. Phys. C39 (1988) 69. 

[12] T. Csorgo and B. Lorstad, Phys. Rev. C54 (1996) 1396. 

[13] U.A. Wiedemann, P. Scotto and U. Heinz, Phys. Rev. C53 (1996) 918. 

[14] U. Heinz et al., Phys. Lett. B382 (1996) 181; Wu Y.-F. et al., |nucl- 
| th/9607044 , 


[15] B.R. Schlci et al., Phys. Lett. B293 (1992) 275; J. Bolz et ah, ibid. B300 
(1993) 404; Phys. Rev. D47 (1993) 3860. 

[16] U.A. Wiedemann and U. Heinz, [nucl-th/9610043| , subm. to Phys.Rev. 
Lett.; and [nucl-th/9611031| , subm. to Phys. Rev. C. 

[17] U. Heinz, in Correlations and Clustering Phenomena in Subatomic 
Physics, (M.N. Harakeh et ah, eds.), NATO ASI Series B, Plenum, 1996, 
in press (|nucl-th/96090291) . 

[18] E. Shuryak, Phys. Lett. B44 (1973) 387; Sov. J. Nucl. Phys. 18 (1974) 


667. 


[19] H. Heiselberg and A.P. Vischer, |nucl-th/9609022| . 

[20] NA49 Coll., K. Kadija et ah, in 0; H. Appelshauser, PhD thesis, Univ. 
Frankfurt (1996). 

[21] S. Schonfelder, PhD thesis, MPI fiir Physik, Miinchen (1996). 

[22] F. Karsch, lectures given at Kloster Banz, Oct. 8-10, 1996; and this vol¬ 


ume. 

















